Counting statistics and super-Poissonian noise in a quantum dot: 
Time-resolved measurements of electron transport 
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We present time-resolved measurements of electron transport through a quantum dot. The mea- 
surements were performed using a nearby quantum point contact as a charge detector. The rates 
for tunneling through the two barriers connecting the dot to source and drain contacts could be 
determined individually. In the high bias regime, the method was used to probe excited states of 
the dot. Furthermore, we have detected bunching of electrons, leading to super-Poissonian noise. 
We have used the framework of full counting statistics (FCS) to model the experimental data. The 
existence of super-Poissonian noise suggests a long relaxation time for the involved excited state, 
which could be related to the spin relaxation time. 

PACS numbers: 



I. INTRODUCTION 

Studies of current fluctuations in conductors are of 
great interest because they give information about the 
charge carriers in the system and their mutual interac- 
tions, complementary to that obtained by the measure- 
ment of the average current In recent years, the 
method of full counting statistics (FCS) has brought 
renewed interest to the field. Using FCS, fluctuations are 
studied by counting the number of electrons that pass 
through a conductor within a fixed period of time. Since 
this gives direct access to the distribution function of the 
fluctuations, not only the shot noise but also higher or- 
der correlations can be extracted. The method has so far 
mainly been used as a theoretical tool for calculating the 
shot noise in various mesoscopic systems [3[. 

For electron transport through quantum dots, the noise 
is typically of sub-Poissionan nature. This is due to 
the Coulomb blockade, which enhances the correlation 
between electrons and thereby reduces the noise 
However, when several channels with different coupling 
strengths contribute to the electron transport, interac- 
tions can lead to more complexprocesses and to an en- 
hancement of the noise [sj, 0,0, Moreover, it has been 
predicted that entangled electrons could lead to super- 
Poissonian noise, thus providing a possible way of de- 
tecting entanglement in mesoscopic systems [qI, To^. 

Experimentally, direct observations of FCS by count- 
ing electrons are difficult to achieve. This is because a 
very sensitive, non-invasive, high bandwidth charge de- 
tector is needed in order to be able to resolve individ- 
ual electrons pjj, ll2|, ll3| . Only very recently, measure- 
ments of FCS for single level transport through a quan- 
tum dot (QD) were performed [IJ. A quantum point 



contact (QPC) was used to read out the charge state 
of the nearby QD [15]. Here, we present further time- 
resolved measurements of a QD system. We show meth- 
ods for tuning the QD and for extracting information 
about tunneling rates [16] and about excited states of 
the QD Furthermore, we present measurements in 
a regime where transport is governed by more complex 
processes than tunneling through a single QD level. We 
observe bunching of electrons and super-Poissonian noise. 
In this regime, we show that the theory of FCS [2] can 
be used to model the experimental data and to extract 
intrinsic properties of the mesoscopic system, such as the 
relaxation time between excited states. 



II. EXPERIMENTAL SETUP 
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FIG. 1: (Color online) (a) Quantum dot with integrated 
charge read-out used in the experiment, (b) Current through 
the QPC as a function of time, showing a few electrons tun- 
neling into and out off the dot. The arrows mark the steps 
corresponding to an electron entering the dot. 

The QD used in the experiment is shown in Fig. 
[Tfa). The structure was fabricated using scanning probe 
lithography [18] on a GaAs/Alo.aCao.rAs heterostructure 
with a two-dimensional electron gas (2DEG) 34 nm below 
the surface (electron density 4.5 x 10~^^ m~^, mobility 
25 m^/Vs). The sample consists of a QD [dotted cir- 
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cle in Fig. [IJa)] and a nearby QPC. We estimate from 
the geometry and the characteristic energy scales that 
the dot contains about 30 electrons. The gates Gl and 
G2 were used to tune the tunnel barriers connecting the 
dot to source and drain leads, while the P gate was used 
to tune the conductance of the QPC to a regime where 
the sensitivity to changes in the dot charge is maximal. 
For our setup, the best sensitivity was reached when the 
QPC conductance (Ggpc) was tuned below the first con- 
ductance plateau, with Gqpc ~ 0.25e^//i. The conduc- 
tance was measured by applying a voltage over the QPC 
(^ias = 500 /iV) and monitoring the current. Since 
changing the voltages on gates Gl and G2 also affects 
the QPC sensitivity, a compensation voltage had to be 
applied to the P-gate in order to keep the QPC in the 
region of maximum sensitivity whenever the other gates 
were changed. All measurements where performed in a 
dilution refrigerator with a base temperature of 60 mK. 

When an electron tunnels onto the dot, the conduc- 
tance through the QPC is reduced due to the electro- 
static coupling between the dot and the QPC. A typical 
time trace of the QPC current is plotted in Fig. [ijb), 
showing switching between two levels. The low levels 
correspond to the configuration where the dot contains 
one extra electron, while Tin and Tout specify the time it 
takes for an electron to tunnel into and out of the dot, re- 
spectively. The length of each time trace presented here 
is 0.5 s. 

The bandwidth of the QPC circuit is A/ = 30 kHz, 
which limits the current we can measure by counting elec- 
trons to / < eAf ~ 5 fA. The bandwidth is similar to 
what was achieved in measurements on a split-gate de- 
fined dot In our setup, the bandwidth is not limited 
by low signal-to- noise ratio (S/N), but by the low-pass 
filter formed by the cable capacitance and the feedback 
resistor of the IV-converter. From the trace shown in Fig. 
[TJb), we extract S/N ~ 15. Assuming a flat noise spec- 
trum, we estimate that S/N would allow us to increase 
the bandwidth by a factor of ten and still get a detectable 
signal. One possible reason for the high sensitivity of 
our detector compared to split-gate defined structures is 
that there are no metallic gates on the surface that shield 
the electrostatic coupling between the QD and the QPC. 
However, the sensitivity is also strongly dependent on 
the exact shape of the confinement potential within the 
QPC and on how susceptible this potential is to changes 
in the electrostatical environment. This may vary a lot 
from sample to sample. For the structure used in this 
measurement, the QPC showed broad resonances in ad- 
dition to the standard plateau features. By operating the 
QPC at the flank of a resonance in the step below the 
first plateau, we were able to find a regime with good 
sensitivity. For typical current levels in the QPC (nA), 
electrons pass the QPC at a rate which is many orders of 
magnitude higher than the rate for electrons passing the 
quantum dot (aA). We conclude that back action of the 
QPC via its shot noise can be neglected for the analysis 
of the counting statistics [2Q] . 



III. THERMAL NOISE WITH ONE LEAD 
CONNECTED TO THE DOT 

In the following, we are interested in the number of 
electrons visiting the dot during a given time interval. 
We call each visit one event and use the symbol r^; to 
denote the number of events occurring per second. In 
the low-bias Coulomb blockade regime, the dot can only 
hold one excess electron. Before a new one can enter, 
another one has to go out. In this case, we can count the 
events by detecting the electrons as they enter the dot 
[marked by vertical arrows in Fig. [Hb)]. Note that by 
counting events, we do not distinguish between electrons 
passing through the dot and electrons hopping back and 
forth between the dot and a single lead. 

First, we concentrate on the regime where only one 
lead is connected to the dot and the electron motion is 
entirely governed by thermal fluctuations and occupation 
probabilities. For the data shown in Fig. [ijb), the gates 
are tuned such that the tunnel barrier between the dot 
and the drain lead is completely closed, while the source 
lead is weakly coupled to the dot. With only one lead 
open and with a temperature and level broadening much 
lower than the charging energy and single level spacing of 
the QD, only one QD state is available for tunneling. The 
probability for an electron to tunnel into or out of the dot 
during a time interval dt is governed by the relation 

Pin/out(^)dt = ri,/oute"^^"/°"^Mt, (1) 

where Fin and Tout are the effective rates for tunneling 
into and out of the dot. Using similar methods as in Ref. 
[13, [l^ , we have checked that Eq. [1] is fulfilled when 
we are in the single-level regime. In the following, we 
consider the case of a non-degenerate level, as discussed 
in previous work [Tg*]. To unambiguously determine the 
spin configuration of the involved states, measurements 
as a function of magnetic field would be required. Given 
the availability of our experimental data we focus our 
analysis on the extraction of the tunneling rates for a 
single non-degenerate state. We also assume the tunnel 
coupling to be independent of energy within the small 
interval of interest. At the end of this section we dis- 
cuss how the analysis would change if spin degeneracy is 
relevant. 

The relations between the effective rates and the dot- 
lead tunnel coupling F are given by 

T,n = rf{AE/kBT), rout = r(i-/(A^/fcBT)), (2) 

where f{x) is the Fermi distribution function, T is the 
temperature and AE is the energy difference between the 
Fermi level of the lead and the electrochemical potential 
of the dot. 

The tunneling rates can be determined directly from 
the measured time traces. Using Eq. (TJ we find Fin = 
l/('?"in), Pout = 1/i^out)^ with a relative accuracy of 
Y^2 \n2/N (see Appendix [A|) . Here, TV is the total num- 
ber of switches occurring during one trace. The rela- 
tive accuracy is calculated assuming that Eq. [1] is valid. 
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For the trace in Fig. mb), we get Fin = 8.0 kHz, 
Fout = 2.3 kHz and F = 1/{t,^) + l/(rout) = 10.3 kHz, 
with a relative accuracy of 1.7%. It has been shown that 
the finite bandwidth of the detector leads to a system- 
atic under-estimate of the actual rates. For the rates 
given here we have compensated for such errors using 
the methods presented in Ref. ^] , with a detection rate 
of Fdet = 100 kHz. 




FIG. 2: (Color online) (a) Average dot population versus volt- 
age on gate G2. The data was fit to a Fermi distribution 
function with T — 230 mK. (b) Counts of events per second 
for the same data as in (a). The data was fit to Eq. (|4|, 
giving F = 2.63 kHz and T — 230 mK. (c) Events per second 
versus Vgi and Vg2- For low values of Vbi and Vb2, both 
the source lead and the drain lead are pinched off. For high 
voltages, the barriers open up so much that the tunneling oc- 
curs on a timescale faster than the measurement bandwidth, 
(d) Temperature (squares) and tunnel coupling (crosses), ex- 
tracted from data shown within the ellipse in (c). As Vg2 is 
increased, Vgi is decreased, in order to keep the dot at a con- 
stant potential. For low Vg2, tunneling occurs between the 
source lead and the dot, for high Vg2, the electrons tunnel 
between the drain and the dot. For intermediate gate values, 
both leads contribute to the tunneling. The electron temper- 
ature was found to be the same for both leads, within the 
accuracy of the data analysis. 

Since the dot can only hold one extra electron, we can 
determine the Fermi function from the average popula- 
tion of excess electrons on the dot 

f{AE/kBT) = (nexcess) = (rout)/((rin) + (Tout)). (3) 

The Fermi function can also be found by counting the 
average number of events occuring per second, te- As- 
suming sequential tunneling and using Eq. ([2]), we find 
for the case with one lead open 



rB = l/((Tin) + {Tout» = r/(l-/). 



(4) 



to change the electrochemical potential of the dot. The 
accuracy obtained when determining the Fermi function 
is A/ = /(I — f)^2/N (see Appendix [B]) . giving error 
bars smaller than the markers used in the figures. The 
data fits well to the expected relations. By first determin- 
ing the lever arm between gate G2 and the dot from stan- 
dard Coulomb diamond measurements [jH, it was possi- 
ble to extract the electronic temperature (T = 230 mK) 
from the width of the Fermi function. The same tem- 
perature was found by checking the width of standard 
Coulomb blockade current peaks [22], measured when the 
dot was in a more open regime. 

As mentioned earlier, the results shown here are valid 
only for a non-degenerate level. Taking spin degeneracy 
into account will modify the value extracted for the tun- 
neling rate F (see Appendix [A|). but it will not change 
the width of the Fermi distribution. The error analysis 
is performed within the assumption of a non-degenerate 
level. If spin-degeneracy is taken into account, then some 
tunneling rates are changed by a factor of two, while the 
expression for the relative accuracy remains the same 
as derived before. Further experiments are required to 
clearly differentiate between spin-degenerate levels and 
single levels at the Fermi energy. However, the general 
way how our analysis proceeds is not affected by this. 



IV. THERMAL NOISE WITH TWO LEADS 
CONNECTED TO THE DOT 

In order to perform time-resolved measurements of 
electron transport through the dot, the tunnel barriers 
have to be symmetrized so that both give similar tunnel- 
ing rates. The rates must be kept lower than the band- 
width of the setup, but still high enough to give good 
statistics. Figure Wis) shows the number of events per 
second as a function of the two gates Vqi and Vg2' In 
the upper left corner of the figure, Vgi is high and Vg2 
is low, corresponding to the case where the source lead 
is open and the drain lead is closed. In the bottom right 
corner, the opposite is true. For the region in between, 
marked by the ellipse in Fig. Efc), the data indicates that 
both leads are weakly coupled to the dot. 

The measurement method does not enable us to dis- 
tinguish whether an electron that tunnels into the dot 
arrives from the left or from the right lead. Therefore, 
when both leads are connected to the dot, the rates in 
Eq. ([2]) must be adjusted to contain one part for the left 
lead and one part for the right lead. 
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= rL(i-/L) + rfl(i-/fl). (5) 



In Fig. [2](a) and (b) we plot the average population and 
the number of events per second as the gate G2 was used 



Here, /l and /r are the Fermi distribution functions of 
the left and the right lead, respectively. Using Eq. ([5]), 
we calculate the rate of events for the case when both 
leads are kept open, 

[Fl/l + TRfR][rL{l - h) + Tnil - fn) 



(6) 



With no bias applied to the dot, the two distributions 
functions /l and fji are identical except for a possible dif- 
ference in electronic temperature in the two leads. How- 
ever, assuming Tl = Tr = T, we have fL^fn^f^ and 
Eq. ^ simplifies to r^; = (Fl + Tr) /(I - /). Fitting 
this expression to curves similar to that shown in Fig. 
[2]^b), we extract the temperature and combined tunnel- 
ing rate F^ + F^^ from the data within the ellipse of Fig. 
Efc). The result is presented in Fig. [2jd). The rates and 
the temperature shown in the graph are due to the com- 
bined tunneling through both leads. Still, for low Vg2 
(high Fgi), the drain lead is pinched off and tunneling 
occurs mainly between the source lead and the dot. For 
high Vg2 (low Fgi), the source is pinched off and the 
tunneling is dominated by electrons going between the 
drain and the dot. The fact that the electronic temper- 
atures extracted from both regimes turn out to be the 
same within the accuracy of the analysis (T=230 mK) 
justifies the assumption that Tl = Tr. 



V. SHOT NOISE AT FINITE BIAS 

Now we apply a finite voltage bias between source and 
drain leads and measure electron transport through the 
dot. Figure ma) shows the Coulomb blockade diamonds 
measured by counting events. In this measurement, the 
gate Gl was used as a plunger gate to control the dot elec- 
trochemical potential. However, the gate also strongly 
affects the source tunnel barrier. For low Gl voltages, 
the source lead is closed, giving strong charge fluctua- 
tions only when the drain lead is in resonance with the 
dot [see case I in Fig. [3l^a,b)]. 

At higher gate voltages, the source lead opens up and 
a current can flow through the dot. In point II of Fig. 
[Sfa), the dot electrochemical potential fin lies within the 
bias window but far away from the thermal broadening 
of the Fermi distribution in the leads. The condition can 
be expressed as 



\±eV/2-fin\ -^ksT, 



(J) 



where the case refers to the source contact and the 
"-" case refers to the drain. Whenever Eq. ([7j) is fulfilled, 
electrons can only enter the dot from the source lead and 
only leave through the drain. In this regime, we measure 
the current through the dot by counting events. This 
opens the possibility to use the QD as a very precise cur- 
rent meter for measuring sub-fA currents [13]. Since the 
electrons are detected one by one, the noise and higher 
order correlations of the current can also be experimen- 
tally investigated [13]. In this regime we measure the 
shot noise of the system, which arises because of the dis- 
creteness of the charge carriers Q. This is in contrast 
to the results shown in the previous sections, where the 
fluctuations were due to thermal effects. 

When the bias exceeds the dot charging energy, Ec ^ 
2.1 meV, and the electrochemical potentials of the (n) 
and the (n + 1) states are within the bias window [see 
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FIG. 3: (Color online) (a) Coulomb diamonds, measured by 
counting events per second. For low values of Vbi, the source 
lead is pinched off and tunneling can only occur between the 
dot and the drain lead. As Vgi increases, the source lead 
opens up and a current can flow through the dot. (b) Dia- 
grams depicting the energy levels of the dot at points I, II and 
ni. In case III, the bias is higher than the charging energy of 
the dot, meaning that the dot can contain 0, 1 or 2 excess 
electrons, (c) Time trace taken at point III. The three possi- 
ble dot populations (n, n + 1 or n + 2 electrons) are clearly 
resolvable. 



case HI of Fig [3l^a,b)], transport processes are allowed 
where the dot may contain 0, 1 or 2 excess electrons. A 
time trace measured at point III of Fig. [31(a) is shown 
in Fig. [3](c). The high sensitivity of the QPC charge 
detector allows us to measure switching between three 
different levels, corresponding to (n), (n + 1) and (n + 2) 
electrons on the dot. This distinction is not possible in a 
standard current measurement. 

With the condition given by Eq. ([7j) fulfilled, we know 
that for positive bias voltage, electrons always enter the 
dot through the source contact and leave the dot through 
the drain contact. In this case, we have 



l/(rin), Fz) = Fo 



1/(to 



(8) 



Equation ([8]) can then be used to determine the tunnel- 
ing rates of an individual state, but only if there are no 
excited states available within the bias window. If there 
are excited states available, Eq. ([8|) will still be valid, 
however, the calculated F^ and Vd will not be the tun- 
neling rates of a single state but rather the sum of rates 
from all states contributing to the tunneling process. A 
further complication with excited states is that there may 
be equilibrium charge fluctuations between the lead and 
the excited state, thereby removing the unidirectionality 
of the electron motion. However, if the relaxation rate of 
the excited state into the ground state is orders of mag- 
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nitude faster than the tunnehng out rate, the electron in 
the excited state wih have time to relax to the ground 
state before equilibrium fluctuations can take place. 

The separate rates Fin and Tout for a close-up of the 
upper-left region of Fig. [Sja) are plotted in Fig. [H^a) 
and (b). It is important to note that the requirement of 
Eq. ([7j) is met only for the region along and above the 
dashed lines in the figures. At the lower left end of the 
dashed lines, the energy levels of the dot are aligned as 
shown in Fig. W^c). Going diagonally upward along the 
lines corresponds to raising the Fermi level of the source 
lead, while keeping the energy difference between the dot 
and the drain lead fixed. 




FIG. 4: (Golor online) (a) and (b): Blow-up of the upper left 
region of Fig. EJa), showing the rates for electrons tunneling 
in (a) and out (b) of the dot, respectively. The solid lines 
mark the positions where the source lead is lining up with 
the electrochemical potential of the dot's ground state. The 
dashed lines mark the lower edge of the region where condition 
of Eq. (|7|) in the text is fulfilled. The color scales are different 
for the two figures, the rate for tunneling out is around 10 
times faster than tunneling in. (c) Diagram depicting the 
energy levels along the dashed lines in (a) and (b). As the 
source lead is raised [corresponds to going upward along the 
dashed lines in (a) and (b)], excited states become available 
for tunneling, (d) Tunneling rate for electrons entering the 
dot, measured along the dashed line in (a). With increased 
gate voltage, excited states become available for transport, 
giving higher tunneling rates. 

Starting at low bias and low voltage on the gate Vgi, 
the dot is in the Coulomb blockade regime, and no tun- 
neling is possible. Following the dashed line upwards, 
the dot ground state becomes available for tunneling at 
^ias = 0.3 mV. The transitions is marked by the solid 
lines in Fig. [3Ka,b). At these low gate voltages, the source 
tunnel barrier is almost completely pinched off, meaning 
that the rate for electrons entering the dot is still low 
[Fig. [41(a)]. Even so, some electrons do enter the dot, 
as can be seen from the few points of measurements of 



rates for electrons tunneling out of the dot within the 
corresponding region of Fig. [11(b). 

We now concentrate on the tunneling-in rate in Fig. 
[11(a). As the source level is further raised, excited states 
become available for transport. The first excited state (at 
"^bias = 0.85 mV along the dashed line) is more strongly 
coupled to the lead than the ground state, giving a tun- 
neling rate of ~ 70 Hz for electrons entering the dot. 
The large difference in the tunneling-in rate between the 
ground and the excited state can be understood if the 
wavefunctions of the ground and excited state have dif- 
ferent spatial distributions. If the overlap with the lead 
wavefunction is larger for the excited state, the tunneling 
rate will also be larger. Similar differences in tunneling 
rates have been found between the singlet and triplet 
states in a two-electron dot [23|, [IJ . 

By further raising the source level, tunneling can 
also occur through a second excited state. The mea- 
sured tunneling-in rate will now be the sum of the rates 
from both excited states; by subtracting the contribu- 
tion from the first state, the tunneling-in rate for the 
second state can be determined. Using this method, we 
can resolve three excited states, with excitations energies 
ei = 0.55 meV, 82 = 1.0 meV, 83 = 1.3 meV and with 
tunneling rates Fi = 70 Hz, F2 = 190 Hz, F3 = 190 Hz. 
The excited states are clearly seen in Fig. [IJd), which is 
a cut along the dashed diagonal line in Fig. [l^a). 

Focusing now on the rates for electrons tunneling out 
of the dot [Fig. [1Kb)], there is a noisy region where the 
ground state but no excited states are within the bias 
window (0.3 < Vhias < 0.85 mV along the dashed line). 
In this regime, few electrons will enter the dot, meaning 
that the statistics needed for measuring the rate of elec- 
trons leaving the dot is not sufficient. However, for bias 
voltages higher than the first excited state, the tunneling- 
out rate remains constant along the dashed line. This is 
in contrast to the steps seen in the tunneling-in rates, 
indicating that the rate for tunneling out of the QD does 
not depend on the state used for tunneling into the QD. 
Since the individual excited states are expected to have 
different rates also for tunneling out of the dot, the data 
is consistent with the interpretation that an electron en- 
tering the dot into an excited state will always have time 
to relax to the ground state before it tunnels out. The 
rate for tunneling out is ~6 kHz, giving an upper bound 
for the relaxation time of ~170 /is. 

The main relaxation mechanism in quantum dots is 
thought to be electron-phonon scattering [25]. Measure- 
ments on few-electron vertical quantum dots have shown 
relaxation times of 10 ns [2^. Recent numerical investi- 
gations have shown that the electron-electron interaction 
in multi-electron dots can lead to reduced relaxation rates 
[27|]. Still, the relaxation rate is expected to be consider- 
ably faster than the upper limit we give here. 
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VI. BUNCHING OF ELECTRONS 

So far, we have analyzed data where the tunnehng 
events can be weh explained by a rate equation approach 
with one rate for electrons tunneling into and another 
rate for electrons leaving the dot. For the trace shown 
in Fig. [51(a), the behavior is distinctly different. The 
electrons come in bunches; there are intervals where tun- 
neling occurs on a fast timescale (>10 kHz), in-between 
these intervals there are long periods of time (> 1 ms) 
without any tunneling. The data was taken with a bias 
applied so that the Fermi level of the source lead is lining 
up with the electrochemical potential of the dot, while 
the drain lead is far below the electrochemical potential 
of the dot, thus prohibiting electrons from entering the 
dot from the drain lead. The voltage on gate Vqi was set 
to 34 mV, which is outside the range of the Coulomb di- 
amonds presented in Fig. [3](a). Since the QPC current is 
at the high level during the intervals without tunneling, 
the dot contains one electron less when the fast tunneling 
is blocked. 
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FIG. 5: (Color online) (a) Time trace of the QPC current 
showing bunching of electrons, (b) Dot states included in the 
model used to describe the bunching of electrons. The red 
circles correspond to electron occupation. State Sa is the n- 
electron ground state, state Sb is an excited n-electron state 
and state Sn+i is the ground state when the dot contains 
(n+1) electrons, (c) Energy diagram for the model. The two 
dot transitions are both within the thermal broadening of the 
lead. Electrons enter the dot from the left lead and may leave 
through either the left or the right lead, (d) Possible tran- 
sitions between the different states of the model. The rates 
Pi^? Pki refer to electrons entering the dot, thus taking the 
dot from state Sa/b to state Sn+i- The rates P^ut? Pout de- 
scribe electrons leaving the dot, giving transitions from state 
5'n+i to Sa/b- Wab and Wba are the direct transition rates 
between states Sa and Sb- Finally, the rates P^ght? Pdght 
refer to electrons leaving the dot through the right lead. 



In order to explain the two different timescales, we as- 
sume the validity of a model where there are two almost 



energy-degenerate dot states within the thermal broad- 
ening of the distribution in the source lead. Because of 
Coulomb blockade, the dot may hold one or zero excess 
electrons. The model includes three possible dot states, 
shown in Fig. [5](b). State Sa is the n-electron ground 
state, state Sb is an excited n-electron state and state 
Sn+i is the ground state when the dot contains (n + 1) 
electrons. Transitions between the Sa/Sb states and the 
6'n+i state occur whenever an electron tunnels into or 
out of the dot. 

The tunnel coupling between the dot and the lead 
is given by the overlap of the dot and lead electronic 
wavefunctions. Since the wavefunctions corresponding 
to the two states Sa and Sb may have different spa- 
tial distributions, the coupling strength P^ of the tran- 
sition Sa ^ Sn-\-i can vary from the coupling P^ of the 
Sb ^ 5'n+i transition. The energy levels of the dot and 
the leads for the configuration where we measure bunch- 
ing of electrons are shown in Fig. El^c), while the possible 
transitions of the model are depicted in Fig. [5](d). 

Starting with one excess electron on the dot [state 5'n+i 
in Fig. EJd)], at some point an electron will tunnel out, 
leaving the dot in either state Sa or state Sb- Assuming 
Tb ^ Pa, it is most likely that the dot will end up in the 
excited state Sb- If the tunneling rate Pb is faster than 
the relaxation process Sb =^ 5^, an electron from the 
lead will have time to tunnel onto the dot again and take 
the dot back to the initial Sn-\-i state. The whole process 
can then be repeated, leading to the fast tunneling in 
Fig. Ha). 

However, at some point the dot will end up in state 
Sa^ either through an electron leaving the dot via the 
Ta transition, or through relaxation of the Sb state. To 
get out of state Sa^ there must be either a direct transi- 
tion back to state 5^, or an electron tunneling into the 
dot through the Sa =^ 5'n+i transition. With P^ ^ Pa 
and assuming Pb ^ Wba^ both processes are slow com- 
pared to the tunneling between the lead and state Sb- 
This mechanism will block the fast tunneling and pro- 
duce the intervals without switching events seen in Fig. 
[5l^a). Similar arguments can be used to show that the 
blocking mechanism will be possible also if P^ <C Pa- 

From the above reasoning, we see that the fast 
timescale is set by the fast tunneling state, while the slow 
timescale is determined either by the relaxation process 
5b ^ 5a or by the slow tunneling rate, depending on 
which process is the fastest. Either way, it is crucial that 
the relaxation rate is slower than the fast tunneling rate 
(in our case Wab <^ Pb ~ 20 kHz). We speculate that 
the slow relaxation rate may be due to different spin con- 
figurations of the two states. For a few-electron QD, spin 
relaxation times of Ti > 1 ms have been reported [23, [23| • 

To make quantitative comparisons between the model 
and the data, we use the framework of full counting 
statistics (PCS) to investigate how the dot charge fluctu- 
ations change as the source lead is swept over a Coulomb 
resonance. Theoretical investigations of multi- level quan- 
tum dots have lead to predictions of electron bunching 
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and siiper-Poissonian noise 0. Following the lines of 
Refs. [3, [2^, we first write the master equation for the 
system, 



(9) 









\ = mITb I 


\ Pn+1 J 


V Pn+l J 



with M = 



-Tf, - Wba 
Wba 



Wab (r 
-rg-WAB {r. 



out 
B 
out 



right 
^ right 
-Tout 



) 

) *e^^ 



Here To 



(r^ut " 



(10) 

- rfut + ^ight + r^ght) and PA, PB and 
Pn-\-i are occupation probabilities for states Sa and 
and 5'n+i, respectively. The effective tunneling rates are 
determined by multiplying the tunnel coupling constants 
for each state with the Fermi distribution of the electrons 
in the lead, 



.A/B 
in/out 



f[T{eV - ija/b)]Fa/b- 



(11) 



The tunneling rates F^gj^^ and F^gj^^ are included to ac- 
count for the possibility for electrons to leave through 
the right barrier. The Fermi level of the right lead is far 
below the electrochemical potential of the dot, so that 
the states in the right lead can be assumed to be unoc- 
cupied. Finally, Wab and Wba are the direct transition 
rates between states Sa and Sb- These rates obey de- 
tailed balance. 



Wab /Wba = exp [{/ha - I^B)/kBT] . 



(12) 



The phenomenological relaxation rate between the two 
states is given as 1/Ti = Wab + Wba- 

In Eq. (p!Q|), we introduce charge counting by multi- 
plying all entries of M involving an electron leaving the 
dot with the counting factor exp(ix) [29]. We do not 
distinguish whether the electron leaves the dot through 
the left or the right lead. In this way we obtain the 
counting statistics pt^iN)^ which is the probability for 
counting N events within the time span to. The distri- 
bution describes fluctuations of charge on the dot, which 
is exactly what is measured by the QPC detector in the 
experiment. We stress that this distribution is equal to 
the distribution of current fluctuations only if it can be 
safely assumed that the electron motion is unidirectional. 
This is the case if the condition in Eq. ([7j) is fulfilled, i.e. 
if the tunneling due to thermal fluctuations is suppressed. 
Here, we are in a regime where there is a mixture of tun- 
neling due to the applied bias and tunneling due to equi- 
librium fluctuations. But since the model defined in Eq. 
(p!Q|) is valid regardless of the direction of the electron mo- 
tion, it can still be used for analyzing the experimental 
data. 

Using the method of Ref. [29], we calculate the lowest 
eigenvalue Aq (x) of M and use it to obtain the cumulant 
generating function (CGF) for p^q(A/'), 



The CGF can then be used to obtain the cumulants of 
any order using the relation Cn = —{—id^)'^S{x)\x^o- 
In order to compare the theory with the experiment we 
extract the first three cumulants of pt^ {N) from the ex- 
perimental data. Since we want to compare the data 
with the predictions given by the CGF of the model, we 
choose to calculate the cumulants instead of the central 
moments, as it was done in a previous work 0]. The 
first cumulant (Ci) is identical to (A^), the mean of the 
distribution, while the second and third cumulants (C2, 
Cs) coincide with the second and third central moments 
[(TV^) - (7V)2 and {N - (N))^], giving the variance and 
the asymmetry of the distribution. 

The cumulants were found by taking a trace of length 
T = 0.5 s and splitting it into m = T/to independent 
traces. By counting the number of electrons N leaving 
the dot in each trace and repeating the procedure for all 
m sub-traces, the distribution function ptQ{N) could be 
experimentally determined. The experimental cumulants 
were then calculated directly from the measured distri- 
bution function [sO]- The time to was chosen such that 
(N) « 3. 
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FIG. 6: (Color online) (a) First, second and third cumulants 
of the distribution of charge fluctuations. The markers show 
values extracted from the experimental data, while the solid 
lines are calculated from the model given in the text. Fit- 
ting parameters are: Fa — 1.6 kHz, Vb — 20.5 kHz, F^g^t = 
4.6 kHz, F^ght = 310 Hz, Ti = 8 ms and ijla - ^jlb = 13 /xeV. 
The electronic temperature was 400 mK. (b) Normalized cu- 
mulants C3/C1 and C2/C1 versus bias voltage. The noise is 
clearly super- Poissonian in the central region of the graph, 
(c) Calculated maximal value of C3/C1 as a function of the 
relaxation time between the two states. The values are calcu- 
lated by varying the relaxation time while keeping the other 
parameters to the values given by the fit shown in (a). The 
maximum value C3 / Ci extracted from the experimental data 
is 15.9. 

Figure [61(a) shows the first three cumulants versus 
voltage applied to the source lead. The points corre- 
spond to experimental data, while the solid lines show 
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the cumulants calculated from the CGF of our model, 
with fitting parameters Ta = 1-6 kHz, Tb = 20.5 kHz, 
Tright = 4.6 kHz, r^ght = 310 Hz, Ti = 8 ms and 
IJ^A — l-^B = 13 /ieV. The electronic temperature in this 
measurement was 400 mK. The figure shows good agree- 
ment between the model and the experimental data. 

Figure [61(b) shows the normalized cumulants C2/C1 
and C3/C1 for the experimental data; we notice that 
both the second and the third cumulants vastly exceed 
the first cumulant when the Fermi level of the source 
lead is aligned with the electrochemical potential of the 
dot (T^Dias = 1.3 mV). The noise is of super-Poissonian 
nature, as expected from the bunching behavior of the 
electrons. 

When the bias voltage is further increased (Vbias > 
1.5 mV), the source lead is no longer in resonance with 
the electrochemical potential of the dot and the equi- 
librium fiuctuations between the source and the dot are 
suppressed. In this regime, the measured charge fiuc- 
tuations are due to a current fiowing through the dot. 
Electrons enter the dot from the source lead and leave 
the dot through the drain lead. The blocking mechanism 
is no longer effective and the transport process will pre- 
dominantly take place through state Sa^ since the tun- 
nel coupling to the drain lead is stronger for this state 
(^right ^ r^gj^^). The transport through the dot can es- 
sentially be described by a rate equation, with one rate 
for electrons entering and another rate for electrons leav- 
ing the dot. For such systems, it has been shown that the 
Coulomb blockade will lead to an increase in correlation 
between the tunneling electrons compared to a single- 
barrier structure, giving sub-Poissonian noise [4, 14]. The 
effect is seen for Fbias > 1-5 mV in Fig [61(b); both the 
second and third cumulants are reduced compared to the 
first cumulant. 

The value of Ti = 8 ms obtained from fitting the ex- 
perimental data is of the same order of magnitude as pre- 
viously reported values for the spin relaxation time Ti. 
We stress that the bunching of electrons and the super- 
Poissonian noise can only exist if the relaxation time is 
at least as long as the inverse tunneling time. This is 
demonstrated in Fig. [He), which shows the maximum 
value obtained for the ratio C3/C1 calculated for differ- 
ent Ti while keeping the rest of the fitting parameters at 
the values given in the caption of Fig. [6l 



VII. CONCLUSION 

In this work, we have shown that a quantum point con- 
tact can be used for measuring time-resolved transport 
through a weakly coupled quantum dot. The detection 
method allows us to determine the tunneling rates for 
electrons entering and leaving the dot separately. Com- 
paring the different tunneling rates, information about 
the excited states and their relaxation times could be 
extracted. We have shown that the framework of full 
counting statistics together with time-resolved measure- 



ment techniques can be used as a tool for extracting infor- 
mation about electron transport properties of solid state 
systems. 

APPENDIX A: STATISTICS OF TUNNELING 
RATES 

In the single-level regime, the process of an electron 
tunneling into or out of the dot is described by the rate 
equation 

Pin/out (^) = -rin/out Pin/out (^)- (Al) 

Here, Pin/out(^) is the probability density for an electron 
to tunnel into or out of the dot at a time t after a com- 
plementary event. Since the expressions for electrons en- 
tering and leaving the dot are the same, we drop the 
subscripts {in /out) and use the notations p(t) and F to 
describe either one of the two processes. Solving the dif- 
ferential equation and normalizing the resulting distribu- 
tion gives 

p(t)di = Fe-^Mt. (A2) 

Equation IA2I is valid assuming non-degenerate states. 
For the case of spin degeneracy, the rate for tunneling 
into the dot should be multiplied with a factor of two if 
both spin states are initially empty. 

In the experiment, we measure a time trace contain- 
ing a sequence of tunneling times r/c, k = 1, 2, 3, . . . To 
estimate F and its relative accuracy from such a se- 
quence, we need to calculate the probability distribu- 
tion for extracting a certain value F, given a fixed se- 
quence of tunneling times. We start by dividing the 
time axis into bins of width Ar and number them with 
z = 0,l,2,...A tunneling event Tk will be counted in bin 
i if zAr < T/e < (z + 1) Ar. Using Eq. (|A2p and assuming 
Ar <C 1/F, we find that the probability to get a count in 
bin i for a given value of F is equal to 

p(z|F) = FAre-^^"\ (A3) 

A certain sequence {in} is realized with probability 

N 
n=l 

= (rAr)^e"^^^^-o^^^ 

= (FAr)^e-r^^^^*>. (A4) 

Here, rii is the number of times an event falls into bin i, 
Xli^o = ^ is the total number of events in the trace 
and {i) = Xli^o i^ average of z. A certain set of 
bin occupations {rii} can be achieved with many different 
{zn}-series, namely NI/YI'^q nj. Assuming that they all 
occur with the same probability p({zn}|F), we find 

p{{n,}\T) = .^^(rAr)^e-^^-^«. (A5) 
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This is our sampling distribution. For an estimate of F 
we use Bayes theorem 



p(r|{nO)=p(r) 



Pi{ni}\r) 
Pi{ni}) 



(A6) 



Because we have no information on the prior probabihties 
p{r) and p({n^}), the principle of indifference requires 
them to be constants, giving 



p(r|K}) = c(rAT)^e 



(A7) 



where C is constant. Normalization p{T\{ni})dT = 1 
leads to 



Pmni}) 



m 



m 



(T)(r(T))^e-^r(r>^ 



The most likely value of F is therefore F* = 1/ (r) . The 
relative accuracy of this estimate is given by the width 
of the distribution. Setting x = F(r) and evaluating the 
width at half maximum gives 



\n{x) = x-1- ^ln(2). 



(A9) 



For large N we can expand \n{x) in a Taylor series around 
X = 1. Keeping only the first two terms, it follows 

\(x-lf = lln{2) 



Thus the relative accuracy is AF/F = \n{2)/N. 



APPENDIX B: FERMI-DISTRIBUTION 

With just one lead connected to the QD, we adopt 
the model Fin = Dn+iTfiAE/kBT); Font = DnTll - 
f{AE/kBT)]. Here, / is the Fermi distribution of the 
lead, F is the tunnel coupling between the QD and the 
lead, Dn is the spin degeneracy of the state with n elec- 
trons on the QD and AE is the energy difference be- 
tween the electrochemical potential of the QD and the 
Fermi level of the lead. The tunnel coupling is assumed 
to be independent of energy within the small interval of 
interest. Combining the equations gives 



n+l 



(Bl) 



In the following we will focus on the special case of non- 
degenerate states, with Dn = ^n+i = 1- To deter- 
mine the Fermi distribution from a sequence of tunneling 
events, we follow the lines of Appendix [Al and divide the 
time axis into bins of width Ar. The bins are labeled 
with i = 0,1,2,..., and an event will be counted in bin 
i if iAr < Tin/out < + l)Ar. We collect the tunnel- 
ing times in two sets of bins, one for electrons tunneling 
into the dot {{ui})^ and one for electrons leaving the dot 
{{rrii}). Assuming that the events corresponding to tun- 
neling in and tunneling out are uncorrelated and using 
the results of Eq. (|A5p . we get 



p{{ni; mi}|Fin; Fout) = p{{ni}\Tir,)p{{mi}\T out) 
oc (FinFoutAr2)^e-^(^^^^"^^>+^°-^"°->). (B2) 



(A8) 

The principle of indifference gives 



p(Fin;Fout|{ni;mi}) ex 

(FinFoutAr2)^e-^(r-^^->+^°-^^°->). (B3) 

Inserting the result of Eq. (|Bip leads to 



p{f;Tont\{ni;mi}) ex 



Marginalization gives 



p{f\{ni;mi}) = / Foutll^i; ^iD^Fout cx 

^0 

(^)'^(2iV+l)! 



^(T^{Tin) + (Tout))A^^ 

The distribution has a maximum for 



21VTT- (B5) 



/ _ (Tout) N 
l-f (Tin) N+l' 

In the limit of A/" ^ oo, the most likely / is given by 



/ 



(Tout) 



(Tin) + (Tout) ' 



(B6) 



We find the width of p{f\{ni;mi}) by approximating it 
around its maximum by a Gaussian. The result is 



a = /(l-/)v/27iV. 



(B7) 
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